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In this paper, we report our results on charged hadron multiplicities of heavy quark 
initiated jets produced in high energy collisions. After implementing the so-called dead 
cone effect in QCD evolution equations, we find that the average multiplicity decreases 
significantly as compared to the massless case. Finally, we discuss on the transverse 
momentum distribution of light quark initiated jets and emphasize on the comparison 
between our predictions and CDF data. 

1. Introduction 

High-p t jets can be initiated either in a short-distance interaction among partons 
in high energy collisions such as pp, pp, in the DIS e p, the e + e _ annihilation and 
via electroweak (or new physics) processes. 

In this framework, we compute the average (charged) multiplicity of a jet initi- 
ated by a heavy quark. For this purpose, we extend the modified leading logarithmic 
approximation (MLLA) evolution equations ^ to the case where the jet is initiated 
by a heavy (charm, bottom) quark. The average multiplicity of light quark jets 
produced in high energy collisions can be written as N(Y) cx exp | j Y j(y)dy^ (Y 
and y are defined in section [2]) , where 7 ~ 70 + A7 is the anomalous dimension 
that accounts for soft and collinear gluons in the double logarithmic approximation 
(DLA) 70 ~ \fal, in addition to hard collinear gluons A7 ~ a s or single logarithms 
(SLs), which better account for energy conservation and the running of the coupling 
constant a s 

The inclusive fc^-distribution of particles inside a jet has been computed at 
MLLA accuracy in the limiting spectrum approximation i. e. assuming an in- 
frared cutoff Qq equal to Aqcd (f° r a review, see also^). MLLA corrections, of 
relative magnitude 0(^/a^) with respect to the leading double logarithmic approx- 
imation (DLA), were shown to be quite substantial Therefore, we also included 
corrections of order 0(a s ), that is next-to- next-to-leading or Next-to-MLLA (NM- 
LLA). 
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Lastly, both observables, the average multiplicity and the inclusive k±- 
distribution are computed unde r th e assumption of local parton hadron duality 
(LPHD) as hadronization model SHI. 



2. Kinematics and variables 

As known from jet calculus for light quarks, the evolution time parameter deter- 
mining the structure of the parton branching of the primary gluon is given by (for 
a review see HI and references therein) 

y = In (^£) , k± = zQ>Q , Q = E6> Q , (1) 

where k± is the transverse momentum of the gluon emitted off the light quark, Q 
is the virtuality of the jet (or jet hardness), E the energy of the leading parton, 
Qo/E < O < Qo is the emission angle of the gluon (O <C 1), Oo the total half 
opening angle of the jet being fixed by experimental requirements. Let us define in 
this context the variable Y a,s y = Y + hi z, y = In The appearance of this 

scale is a cons eque nce of angular ordering (AO) of successive parton branchings in 
QCD cascades 1221. 

The corresponding evolution time parameter for a jet initiated by a heavy quark 
with energy E and mass m appears in a natural way and reads y = In j K \ = 

k\ + z 2 m 2 1^1, which for collinear emissions 8<1 can also be rewritten in the form 

k ± = zQ, Q = E(e 2 + 2 m )\ (2) 

with O > © m (see Fig.[lJ. An additional comment is in order concerning the AO for 
gluons emitted off the heavy quark. In ([2]), O is the emission angle of the primary 
gluon g being emitted off the heavy quark. Now let O' be the emission angle of 
a second gluon g' relative to the primary gluon with energy u' <C us and 0" the 
emission angle relative to the heavy quark; in this case the incoherence condition 
O' 2 < (O 2 + 0^) together with 0" > m (the emission angle of the second gluon 
should still be larger than the dead cone) naturally leads ((2j) to become the proper 
evolution parameter for the gluon subjet (for more details see SI). For m = 0, the 
standard AO (©' < 0) is recovered. Therefore, for a massless quark, the virtuality of 
the jet simply reduces to Q = EQ as given above. The same quantity k± determines 
the scale of the running coupling a s in the gluon emission off the heavy quark. It 
can be related to the anomalous dimension of the process by 

7o 2 (^) = 2N C ^± = — L_ Mn}) = ±-( l ±N c - 2 -n f ) , A = In ■ °° 



A,(y + A)' ^ V I; 47V C V3" L 3"V Kqcd 

(3) 

where rif is the number of active flavours and iV c the number of colours. The 
variation of the effective coupling a s asn/ — > U/ + 1 over the heavy quarks threshold 
has been suggested by next-to-leading (NLO) calculations in the MS scheme HI and 
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is sub-leading in this frame. In this context (3o( n f) wm be evaluated at the total 
number of quarks we consider in our application. The four scales of the process are 
related as follows, 

Q > m > Qo ~ Aqcd, 

where Qq ~ ^qcd corresponds to the limiting spectrum approximation Finally, 
the dead cone phenomenon imposes the following bounds of integration to the 
perturbative regime 

™<z<l-™, m *<tf<E\<dl + QlX (4) 

which now account for the phase-space of the heavy quark jet. The last inequality 
states that the minimal transverse momentum of the jet Q = EQ m — m is given by 
the mass of the heavy quark, which enters the game as the natural cut-off parameter 
of the perturbative approach. 



3. QCD evolution equations 

The system of QCD evolution equations for the heavy quark initiated jet, corre- 
sponding to the process described in Fig. Q] is found to read^ 

V - ei (Y-,L m ) / dy^(y)N g (y) - e 2 (f, L m ) 7o 2 (f )iV 9 (y), (5) 



C F dY v ' ' Jy, 



dN„ 



! - I "'" d V1 l{y)NM - A(Y,L m h 2 (Y)N g (Y), (6) 



dY j Ym 



where A(Y, L m ), a(rif), £i(Y, L m ) and £2(^1 L m ) are defined in^. As a consistency 




Fig. 1. Parton splitting in the process Q — > Qg: a dead cone with opening angle © m is schematically 
shown. 



check, upon integration over Y of the DLA term in Eq. ([5]) , the phase space structure 
of the radiated quanta reads 

r Q o P) 2 dP) 2 rl-m/Q j 
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Notice that the lower bound over O 2 in ([7} (Y in ((3J)) can be taken down to "0" 
(Y m = L m in iJSJ)) because the heavy quark mass plays the role ol collinear cut-off 
parameter. 

4. Phenomenological consequences 

Working out the structure ol ([5]) and ([5]), we obtain the rough difference between 
the light and heavy quark jet multiplicities, which yields, 



It can be seen that ([5]) is exponentially increasing because it is dominated by the 
leading DLA energy dependence ol N q . According to (jHJ), the gap arising from the 
dead cone effect should be bigger lor the b than for the c quark at the primary state 
bremsstrahlung radiation off the heavy quark jet. The approximated solution ol the 
evolution equations leads to the rough behaviour ol N q — Nq in ([5]), which is not 
exact in its present form. That is why, in the following, we use the numerical solution 
ol the equations © and © ■ In this study we advocate the role ol mean multiplicities 




Fig. 2. Massless and massive quark jet average multiplicity Nj^ tal including heavy quark flavour 
decays as a function of the jet hardness Q. Bands indicate m c and in the [1, 1.5] and [3,5] 
GeV intervals respectively. 

ol jets as a potentially uselul signature for 6-tagging and new associated physics^ 
when combined with other selection criteria. In Fig. [51 we plot as frmction ol the jet 
hardness, the total average jet multiplicity (N^ 111 ), which accounts for the primary 
state radiation off the heavy quark (N^ 1 ) together with the decay products from the 
final-state flavoured hadrons {N^ c ), in the form N^ tal = Nf' + Nj . Moreover, the 
flavour decays constants N^ c — 2.60±0.15 and N{f c — 5.55±0.09^are independent 
ol the hard process inside the cascade, such that Njf can be added in the whole 
energy range. Alter accounting for the weak decay multiplicities, it turns out that 
the b quark jet multiplicity becomes slightly higher than the c quark jet multiplicity, 
although both remain suppressed because ol the dead cone effect. 




1 — exp 




(8) 
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Fig. 3. CDF preliminary results for the inclusive k , distribution at various hardness Q in com- 
parison to MLLA and NMLLA predictions at the limiting spectrum (Qo = ^QCd)\ the boxes are 
the systematic errors. 



5. Single inclusive k±— distribution of charged hadrons in NMLLA 

Computing the single inclusive k±— distribution requires the definition of the jet 
axis. The starting point of our approach consists in considering the correlation be- 
tween two particles (hi) and (h2) of energies E\ and Ei which form a relative angle 
inside one jet of total opening angle Oo > 6 Weighting over the energy E2 of 
particle (h.2), this relation leads to the correlation between the particle (h=hl) and 
the energy flux, which we identify with the jet axis l ™. Thus, the correlation and 
the relative transverse momentum k± between (hi) and (h2) are replaced by the 
correlation, and transverse momentum of (hi) with respect to the direction of the 
energy flux. In Fig. [3]^, we report measure men ts over a wide range of jet hardness, 
Q = EQq, in pp collisions at y/s = 1.96 

TevUni 

together with the MLLA predictions 
of ^ and the NMLLA calculations, both at the limiting spectrum (A = 0) and taking 
Aqcd = 250 MeV; the experimental distributions suffering from large normaliza- 
tion errors, data and theory are normalized to the same bin, \n(k ± /l GeV) = —0.1. 
The hadronic fcj_-spectrum inside a high energy jet is determined including cor- 
rections of relative magnitude O(^Ja^) with respect to the MLLA, in the limiting 
spectrum approximation. The results in the limiting spectrum approximation are 
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found to be in impressive agreement with measurements by the CDF collaboration, 
unlike what occurs at MLLA, pointing out small overall non-perturbative contribu- 
tions. The MLLA predictions prove reliable in a much smaller k ± interval. At fixed 
jet hardness (and thus Yq ), NMLLA calculations prove accordingly trust-able in 
a much larger x interval. 

6. Conclusions 

Thus, our present work focusing on the differences of the average charged hadron 
multiplicity between jets initiated by gluons, light or heavy quarks could indeed 
represent a helpful auxiliary criterion to tag such heavy flavours from background 
for jet hardness Q > 40 GeV. Notice that we are suggesting as a potential signature 
the a posteriori comparison between average jet multiplicities corresponding to 
different samples of events where other criteria to discriminate heavy from light 
quark initiated jets were first applied. Figj2]plainly demonstrate that the separation 
between light quark jets and heavy quark jets is allowed above a few tens of GeV 
with the foreseen errors of the experimentally measured average multiplicities of 
jets. 

The single inclusive /c^-spectra inside a jet is determined including higher-order 
0(a s ) (i.e. NMLLA) corrections from the Taylor expansion of the MLLA evolu- 
tion equations and beyond the limiting spectrum, A ^ 0. The agreement between 
NMLLA predictions and CDF preliminary data in pp collisions at the Tevatron is 
impressive, indicating very small overall non-perturbative corrections and giving 
further support to LPHD ™. 
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